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I. INTRODUCTION 



During the last few years the interest into theoretical investigations of the quark-gluon 
structure of light mesons, especially the pion, rj and rf mesons, has risen due to the high- 
orecision CLEO results on the electromagnetic M7 transition form factors (FF's) Fm^{Q 2 ) 
]J (with M denoting one of these mesons), as well as because of the observed very large 
branching ratios for the exclusive B — > Krj' and semi-inclusive B — > r]'X s decays 0. 

The data on the 7/7 transition FF were mainly used for extracting information concerning 
the 77' quark component of the meson distribution amplitude (DA). Schemes and methods ap- 
plied to this purpose range from light-cone perturbation-theory calculations - with the quark 
transverse- momentum k± dependence kept in the hard-scattering amplitude Th{x, k±,Q 2 ) 
of the underlying hard subprocess - to the modified hard-scattering approach (mHSA) 
with resummed transverse-momentum effects (giving rise to Sudakov suppression factors) 
and such due to the intrinsic transverse momentum of the meson wave functions [H, 111 , to 
the running coupling (RC) method, employed for the estimation of power-suppressed cor- 
rections to the 7/7 transition FF . In these investigations, two different parameterizations 
were used, one employing the conventional 77 — 77' mixing scheme with one mixing angle 9 P 
0- 0- El for both physical states and decay constants and a second one, which considers two 
mixing angles 6\ and 6% to parameterize the weak decay constants f P (P = 77, 77'; 7 = 1,8) 
of the 77 and //' mesons An important conclusion drawn from these investigations, irre- 
spective of the underlying method, is that the 77'-meson DA must be close to its asymptotic 
form and that the admixture of the first non-asymptotic term should be within the range 
62(7^0) — 0-05 0-15, b 2 being the first Gegenbauer coefficient. The CLEO data on the 7/7 
transition and the two-angles mixing scheme were also used to estimate the allowed range 
of the intrinsic charm content of the 77' meson decay constant fZ [5]. It turned out that the 
value —65 MeV <fZ < 15 MeV does not contradict the CLEO data. 

But apart from the ordinary light quark and charm \r)' c ) components, the 77' meson may 
also contain a two-gluon valence Fock state \gg). Moreover, absent at the normalization point 
/Jo, a gluon component of the 77' meson will appear in the region Q 2 > Up owing to quark- 
gluon mixing and renormalization-group evolution of the 77' meson da HUH Ha EH- This 



can directly contribute to the 7/7 transition at the next-to-leading order due to quark box 
diagrams and also affect the leading-order result through evolution of the quark component 
of the T/'-meson DA. Hence, an effect of the 7/'-meson gluon component on the 7/7 transition 
is only mild and was therefore neglected in most theoretical investigations (E El IE El ■ 

But the contribution of the gluon content of the 77' meson to the two-body non-leptonic 
exclusive and semi-inclusive decay ratios of the B meson may be sizeable. Indeed, such a 
mechanism to account for the observed large branching ratio |2| 

Br(B -> T]' + X s ) = (6.2 ± 1.6 ± 1.3) x 10" 4 , (1.1) 

was proposed in Ref. In this work it was suggested that the dominant fraction of 

the B — > 77' + X s decay rate appears as the result of the transition g* — > grj of a virtual 
gluon g* from the standard model penguin diagram b — > sg* . For the computation of the 
contribution of this mechanism to the Br{B — > 77' + X s ) in Ref. the g*gv( vertex 
function was approximated by the constant H(q 2 ,0,m 2 ,) ~ H(0,0,m 2 ,) ~ 1.8 GeV -1 , the 
latter being extracted from the analysis of the J/ip — > 7/7 decay. Further investigations, 
however, demonstrated that effects of the QCD running coupling a s (q 2 ) [13], as well as 
a momentum dependence of the form factor H(q 2 ,0,m 2 ,), properly taken into account, 
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considerably reduce the contribution to Eq. (jl.lj) of the mechanism under consideration 14 1 . 
In order to cover the gap between theoretical predictions and experimental data in Ref . 15j , 
a gluon fusion (spectator) mechanism was proposed. In accordance with the latter, the 7/ 
meson is produced by the fusion of a gluon from the QCD penguin diagram b — > sg* with 
another one emitted by the light quark inside the B meson. In this mechanism, the form 
factor 1 F v > g * g *(qi 2 , gf , m 2 ,) appears owing to the g*g* — > 7/ transition. 

The same ideas form a basis for the computation of the branching ratios of various two- 
body non-leptonic exclusive decay modes of the B meson 0]. To account for the data 
on the exclusive decay B — > Krj in Ref. [13], a mechanism was proposed, based on the 
assumption of a strong intrinsic charm component of the rj' meson. But a more detailed 
analysis 0] proved that the charm content of the rf is too small (f5 ~ —2 MeV) to explain 
the observed branching ratio BriB —>■ Krf). The CLEO data on the B meson non-leptonic 
decays that stimulated interesting theoretical works 0] remain actual until today [2o| . 

The exclusive processes B — > K^g^ were also analyzed within the QCD factorization 



approach [21| and various contributions to the corresponding branching ratios were estimated 
22j| . In accordance with Ref. [22J , the spectator mechanism and such connected to the gluon 
content of the 77 and 7/ mesons are not a key factor in explaining the pattern of the observed 
experimental data. Instead, a more important role plays here the interference (constructive 
or destructive) of penguin amplitudes and large radiative corrections to them, which bring 
the predicted branching ratios into reasonable agreement with the data. 

The //-meson virtual- (on-shell-) gluon transition form factor F v ' g * g *(q 2 , g§, m 2 ,) is the cen- 
tral ingredient of analyses performed within perturbative QCD (pQCD) and hence requires 
further thorough investigations. Such theoretical investigations are especially imp ortant in 
view of contradictory predictions made for this FF in the literature (23I . 0, |25[ |26| . For 
instance, in Ref. |23( the contribution of the gluon content of the 7/ meson to the form factor 
F v igg*(qf , 0, mi/) was estimated and found to be too small while F^gg* (qf, 0, m 2 ,) itself is 
close to the form 

used in some phenomenological applications HQ. 

On the other hand, in computations 
of the //-meson virtual-gluon vertex function F v i g * g * (qf, q%, m 2 ,), the gluon content of the 
7/ meson was neglected and the asymptotic form of the DA for the quark component was 
employed ji^]. The same vertex function F v r g *g*(qf,q%,m 2 ,) was considered in [25], where 
both the standard HSA j2?| as well as the mHSA were used, and the space-, and time-like 
vertex functions were analyzed. Some errors made in Ref. in the description of the hard- 
scattering amplitude of the subprocess q + q — > g* + g* and in the definition of the evolution 
of the 77' meson DA were corrected in 25}. Unfortunately, also this latter work contains an 
incorrect expression for the gluon component F5 g * * (Q 2 , u, g) of the form factor, since it is 
antisymmetric under the exchange of the asymmetry parameter u <-» —uj in clear conflict 
with Bose symmetry. More recently, in Ref. 26], the space-like form factor F rj / g * g *(Q 2 ,uj) 
was computed within the standard HSA. These authors performed a detailed analysis of 
the normalization of the gluon component of both the //-meson DA and that of the gluon 
projector onto a pseudoscalar meson state. 



In this work we use the notions "vertex function" and "form factor" on the same footing. 
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In the present work, we investigate the massless 7/-meson virtual-gluon space-like tran- 
sition form factor F v > g * g * (Q 2 , u) using the framework of the standard HSA |22j], as well 
as by applying the RC method together with the infrared (IR) renormalon calculus |28|. 
Our results obtained within the standard HSA are in agreement with those of Ref. [26]. 
But our central task here is the calculation and evolution of power-suppressed corrections 
~ 1/Q 2n , n = 1,2... to the transition FF Q 2 F r) i g * g *(Q 2 Because in the production of 
the T)' meson from the B decay the momentum squared of the virtual gluon can vary from 
1 GeV 2 to 25 GeV 2 , the power corrections in this domain of Q 2 are expected to play an 
important role. Note, however, that we elaborate only on a theoretical framework to com- 
pute power corrections for the space-like transition FF. Nevertheless, our technique can be 
generalized to encompass the time-like transition form factor, relevant for B meson decays, 
as well. Work in this direction will be reported elsewhere. 

The RC method enables us to estimate power corrections coming from the end-point 
x, y — > 0, 1 regions (for defmiteness we consider two mesons in an exclusive process) in the 
integrals determining the amplitude for an exclusive process. It is known |27[ that in order 
to calculate an amplitude of some hadron exclusive process, one has to perform integra- 
tions over longitudinal momentum fractions x, y of the involved partons. If one chooses the 
renormalization scale fi 2 R in the hard-scattering amplitude Tjj of the corresponding partonic 
subprocess in such a way as to minimize higher-order corrections and allows the QCD cou- 
pling constant a s (n 2 R ) to run, then one encounters divergences arising from the end-point 
x, y — > 0, 1 regions. The reason is that the scale fi R , as a rule, is equal to the momen- 
tum squared of the hard virtual partons, the latter carrying the strong interactions in the 
subprocess' Feynman diagrams |29( and, in general, depends on x and y. Within the RC 
method this problem is resolved by applying the IR renormalon calculus. It turns out that 
this treatment allows us to evaluate power-behaved corrections to the physical quantity un- 
der investigation. This method was recently used for the computation of power corrections 
to the electromagnetic form factors F^ xiO 2 ) of the pion and the kaon, |3(J HH, and the 
electromagnetic transition FF's Fm 7 (Q 2 ) [a |32j of the light pseudoscalar mesons. Power 
corrections to hadronic processes can also be calculated utilizing the Landau-pole free ex- 
pression for the QCD coupling constant [33[. This analytic approach was used to compute in 
a unifying way power corrections to the electromagnetic pion form factor and to the inclusive 
cross section of the Drell-Yan process 0, |35| . 

The paper is organized as follows. In Sect. II we briefly review 77 — 77' mixing schemes in 
the flavor SUf(3) octet-singlet and the quark-flavor bases, and discuss the evolution of the 
quark and gluon components of the ?/ meson DA with the factorization scale. The important 
question of the generalization of the hard-scattering amplitudes of the rj'g* transition to the 
RC method case is also considered. Section III is devoted to a rather detailed presentation of 
the RC method. In Sect. IV we compute the quark and gluon components of the transition 
FF F v i g * g * (Q 2 , uj) . Section V contains our numerical results. Finally, in Sect. VI we make 
our concluding remarks. Appendix contains an additional information on the r/- meson DA. 

II. QUARK AND GLUON CONTRIBUTIONS TO THE rf - g* TRANSITION 
FORM FACTOR 

In this section we consider the quark-gluon content of the rj' meson, as well as the rj 
DA and give some general expressions for determining the rfg* transition form factor within 
both the standard HSA and the RC method. 
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A. Structure of the rf meson 



The parton Fock state decomposition of the pseudoscalar P — rj, rf mesons, can be 
generically written in the following form 

\P) = \P a ) + \P b ) + \P C ) + \P g ) , (2.1) 

where \P a ) and \Pb) denote the P-meson light-quarks u, d, s and |P C ), \P g ) its charm and 
gluon components, respectively. Fock states with additional gluons and qq quark- ant iquark 
pairs have been omitted for simplicity. 

The P-meson light-quark content \P a ) , |P&) can be described either in the flavor SUf(3) 
octet-singlet or in the quark-flavor basis. In the first scheme, the states \P a ) , |P&) are 
expressed as superpositions of the SUf(3) singlet r/i and octet ri 8 states 

*i _ - _ 

\r]i) = -j= uu + dd+ ss 

|r7s> = -p uu + dd-2ss), (2.2) 



whereas in the quark-flavor basis, the SUf(3) strange i] s and non-strange r\ q states are used, 
i.e., 

* -\ 

\Vq) = ^uu + dd), \r} s ) = ^ s \s) . (2.3) 

In Eqs. (12 .2|) and ()2.3|) denote wave functions of the corresponding parton states. 

As mentioned above, the charm component of the rf meson was estimated jlH |36| to be 
too small to considerably affect the B meson exclusive decays. Therefore, in the present 
investigation we neglect the charm content of the rf meson. The maximal admixture of the 
two-gluon state in the rf meson was estimated to be around 26% of its content New 
results from the KLOE Collaboration [3£| are compatible with the two-gluon contribution 
in the rj meson being below the level of 15%. 

The pure light-quark sector of the r\ — r{ system without charm and gluon admixtures 
can be treated by means of the basic states ()2.2|) or f)2.3j) . In the SUf(3) octet-singlet basis, 
the rj and rj' mesons are expressed as superpositions of the rj 8 , rji states 

\r)) = cos P \r) 8 ) - sin 9 P \r]x) , 

\rj') = sin 9 P \rj 8 ) + cos 9 P \r]i) . (2.4) 

Here, 9 P is the mixing angle of physical states in the octet-singlet scheme. In the quark- 
flavor basis we get the same expressions but with 9 P replaced by the mixing angle <p p of the 
physical states in the new basis rj q , ri s , viz., 

\r)) = cos (f) p \rj q ) - sm(p p \r) s ) , 

\r)') = sin0 p \i] q ) + cos^p \rj s ) . (2.5) 

The rji and 7] 8 states and the mixing angle 9 P in the octet-singlet scheme can be expressed 
in terms of the r) q , r) s states and the mixing angle cj) p in the quark-flavor basis and vice versa. 
At the level of physical-state mixing there is no difference between the two bases (12 .2j) and 
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(|2.3|) . This only appears when one parameterizes the decay constants f P of the P = 77 and 
7]' mesons in terms of 

(0\J*b\p) = ifpP* 

where J l 5 is the axial- vector currents with i — q, s or i — 1, 8. In the quark- flavor basis the 
decay constants follow with great accuracy the pattern of the state mixing 3^, 4^ 



fv = fl COS 0P> fv = ~f* Sil1 

/# = /<? sin P> = /a cos <P P - (2-6) 

The situation with the parameterization of the decay constants in the octet-singlet basis is 
different. In this case, in order to take into account the flavor SUf(3) symmetry breaking 



1(8) 

effects, a two mixing-angles scheme for the decay constants / P was introduced [41 

fr, = fa cos fr, = ~fi sin 01, 

# = / 8 sin0 8 , fi= f 1 cos9 1 . (2.7) 

The mixing angles 9\ and 9% differ from each other and also from the state mixing angle 9 P . 
Nevertheless, in some phenomenological applications the conventional octet-singlet mixing 
scheme, that is the scheme with one mixing angle and which assumes the equality 9\ = 
9 8 = 9 p , is used. In the investigation of physical processes both the octet-singlet and the 
quark-flavor bases may be used. A detailed analysis of the rj — rj' mixing problems and 
further references can be found in Ref. 4^ ■ 



In our present work we choose the quark-flavor basis (J2.3|) and the mixing scheme (J2.5j) . 
In this scheme, the decay constants f q and f s , and the mixing angle (fi p have the following 
values HI El 



f q = (1.07 ± 0.02)/,, f s = (1.34 ± 0.06)/,, P = 39.3° ± 1.0° 

with f n = 0.131 GeV being the pion weak decay constant. The KLOE result for <ft p = 
(41.8!j[ 1 6o ) [38] is slightly shifted towards larger values, but it still does not contradict the 
average value <f) p = 39.3° ± 1.0°. In our numerical computations we shall use the central 
values of the constants shown above. 

The singlet part of the rj' meson DA 2 <f>(x,fi 2 ) was obtained in Refs. 0,0,0] by solving 
the evolution equation and found to depend on both functions (p q (x, /z 2 ) and (f> 9 (x, jj 2 ). These 
functions denote the quark and gluon components of the r/ meson DA, respectively, and 
satisfy the symmetry and antisymmetry conditions, 

^(x,/x 2 )=^(l-x,/x 2 ), ^(x,fi 2 ) = -^(l-x,fi 2 ). (2.8) 

This follows from the symmetry properties of the DA of the two-particle bound state of 
a neutral pseudoscalar meson and is in fact enough to obtain general expressions for the 
7/-meson virtual-gluon transition form factor. The evolution of the functions (f) q (x,/j, 2 ) and 
4> 9 (x,fi 2 ) with the scale fi 2 , as well as their dependence on the constants f q , f s and P will 
be considered in Subsect. II. C. 



2 The octet component of the rf meson DA is irrelevant for our present investigation and will not be 
considered here. In what follows we refer to the singlet part of the rj meson DA as being the v( meson 
DA. 
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B. The rf — g* transition form factor F 7? / g * g * (Q 2 , u) 

The massless rf meson - virtual gluon transition form factor F n i g * g * (Q 2 , uj) 

F,, g * g *{Q 2 ^) = F q g * g *(Q 2 ,uj)+F 9 g * g *(Q 2 ,uj), (2.9) 
can be defined in terms of the invariant amplitude 

M = M q + M 9 , (2.10) 

for the process 

r)'(P)^g*( qi )+g*(q 2 ), (2.11) 

in the following way 

MM = -iF^{Q 2 ,u)8^°e a ;e h ; qip q 2a . (2.12) 

In Eq. (12.12)1 e a , e b v and q\, q 2 are, respectively, the polarization vectors and four-momenta 
of the two gluons. Because we study only the space-like FF, q 2 and q\ obey the constraints 
Q\ = -q\ > 0, Ql = -q\ > and Q\ = -q\ > 0, Q\ = -q\ > 0. The form factor 
F v 'g* g *(Q 2 ,u) depends on the total gluon virtuality Q 2 and the asymmetry parameter u, 
defined by 

Q 2 = Ql + Ql " = ^q^- (2-13) 

The asymmetry parameter u varies in the region — 1 < uj < 1 . The value u — ± 1 corresponds 
to the r/-meson on-shell-gluon transition and the value uj = to the situation when the 
gluons have equal virtualities Q\ = Q\. 

In accordance with the factorization theorems of pQCD at high momentum transfer, the 
FF's F^g} g *(Q 2 ,uj) are given by the expressions 

(Q 2 , uj) = [T q H (x, Q 2 , uj, /4) + T q H (x, Q 2 , uj, /£)] ® <P\x, //|) (2.14) 

and 

F 9 v , g , g , (Q 2 , uj) = [T 9 H (x, Q 2 , uj, i4) - T 9 H (x, Q 2 , uj, fi 2 F )} ® <j?{x, fi 2 F ). (2.15) 

Here, x = 1 — x and ftp represents the factorization scale. In Eqs. ()2.14j) and 1)2.15)1 we have 
used the notation 

T H {x, Q 2 ,uj,[i 2 f ) <g> <f)(x, n 2 F ) = [ T h (x,Q 2 ,uj, i&)<f>(x, [j, F )dx. (2.16) 

Jo 

The sum 

T q H (x, Q 2 , uj, hI) + T q H (x, Q 2 , uj, hI) (2.17) 

and the difference 

T 9 H {x, Q 2 , uj, /4) - T 9 H (x, Q 2 , uj, h%) (2.18) 

represent the hard-scattering amplitudes of the subprocesses q + q — > g* + g* and g + g — ► 
g* + g* , respectively. The Feynman diagrams contributing at the leading order to these 
subprocesses' amplitudes are depicted in Figs. [T] and |U 
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FIG. 1: Leading-order Feynman diagrams contributing to the hard-scattering subprocess q + q 
9* +9*- 





FIG. 2: Feynman diagrams contributing at leading order to the subprocess g + g — » g* + g* . 

In what follows, we omit the subscript H in Eqs. (J2.17j) and (|2.18|) and introduce instead 
the notations 

T? 9 \x,Q\u,&)=T$ 9 \x,Q\u>,&), 



Tf 9 \x,Q 2 ,u,^ F )=T^\x,Q",u J ,^ F ). 

In leading order of pQCD, we get for the hard-scattering amplitudes Tf^ and of 
the massless rj' meson - virtual gluon transition FF the following expressions 



T^(x,Q 2 ,uj,fi 2 R ) 



Air 



«s(/4) 



3Q 2 (1 +u)x + (1 -u)x' 
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(X ' ^ ' ^ = -3^ (l + ^+(l-^ (2 ' 19) 



and 

T?(x,Q 2 ,u,[i 2 R ) 



o 2 . 7TQr s (/ij e ) (1 + 0;)x + (1 — Uj)x 



Q 2 rif uj [(1 + uj)x + (1 — u;)x] ' 



^/ n 2 _2n ™s(/4) (1 + Uj)x + (1 ~ Uj)x 

T|(x,Q ,a;,/x fi ) = — = r— — r — — r-r. 2.20 

Q 2 rif uj [{l + uj)x + (1 — cjJxJ 

In deriving the hard-scattering amplitudes Tj^Jx, Q 2 , uj, [i 2 R ) (they are the object of con- 
tradictory conclusions in the literature [HI, the following projection operator of the 



?/-meson onto the two-gluon state has been used 

x 7 S 



fii/,ab 



Note that in the limit uj = the difference Tf — Xf , being the physical hard-scattering 
amplitude, is singularity-free (see Eq. (|4.25|) below). In the expressions above, a s (/i 2 ) is the 
QCD coupling constant in the two-loop approximation defined by 



4tt 



A,ln(// 2 /A 2 ) 



2/3i lnln^/A 2 ^ 
~W Hfi 2 /A 2 



1 - 



(2.21) 



with p Q and Pi being the one- and two-loop coefficients of the QCD beta function: 



2 19 
p = U-~n f , p 1= 51-—n f . (2.22) 

In Eqs. IjOty . fl2~2Tl. and $LT^i A = 0.3 GeV is the QCD scale parameter and n f is the 
number of active quark flavors. 

The physical quantity F rj i g * g *{Q 2 , uj), represented at sufficiently high Q 2 by the factoriza- 
tion formulas (|2.14|) . (J2.15|) . is independent of the renormalization scheme and the renor- 
malization and factorization scales fi 2 R (~p 2 R ) and fi 2 F . Truncation of the perturbation series 
of F n i g * g * (Q 2 , u) at any finite order causes a residual dependence on the scheme as well as 
on these scales. At the leading order of pQCD, the hard-scattering amplitude does not 
depend on the factorization scale fi 2 F but depends implicitly on the renormalization scale fA R 
through a s (fi 2 R ). 3 An explicit dependence of the function Th on fi 2 R and fi 2 F appears at 0(a s ) 
due to QCD corrections. As the scales fj, R and fi F are independent of each other and can 
be chosen separately, we adopt in this work the natural and widely used "default" choice 
for the factorization scale [i F = Q 2 and omit in what follows a dependence of the hard- 
scattering amplitude on fi 2 F (see, Eqs. (I2.19|) and (12.20(0 . For a more detailed discussion of 
these questions, we refer the interested reader to |43j |. 

In the standard HSA with a "frozen" coupling constant j2?J one sets the renormalization 
scale to be fi R = ~p 2 R = Q 2 , simplifying the calculation of the vertex function F ri i g * g *(Q 2 ,uj) 



Similar arguments hold for the scale fi R . 
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considerably. In this approach the hard-scattering amplitudes T x 2 and hence the quark and 
gluon components of the vertex function F n i g * g * (Q 2 , uS) possess symmetry features pertaining 
to the RC method. Indeed, it is not difficult to see that 

T^ 2) (x,Q 2 ,uj) = T 2 9 (1) (x,Q 2 ,u;), T^ 2) (x, Q 2 , -u) = T^ {1) (x,Q 2 ,uj) (2.23) 

and 

Tf (2) (x, Q 2 , u) = T* {1) (x, Q\ u), Tf (2) (x, Q 2 , -u) = -T| (1) (x, Q 2 , u). (2.24) 
Using Eqs. flZB) and (J22H|), we then find 

F} g * g *{Q 2 , to) = T?{x, Q 2 , u) <g> 9 (x, Q 2 ) + {uj^ -u) = 

= 2T?{x,Q 2 ,u)®<i) q {x,Q 2 ). (2.25) 

In similar manner, using Eqs. ()2.8|) and (|2.24j) . one can demonstrate that the following 
equalities hold 

F° gtgt (Q 2 ,u) = Tf(x,Q 2 ,u) <g> ff (x,Q 2 ) + (a; <- -a;) = 

= 2Tf(x,Q 2 ,^)®</> 9 (;r,Q 2 ). (2.26) 
From Eqs. (|2~25J) and flZgj) it follows that 

F« W (QV) = F^ g ,(g 2 , -a;), F^ W (Q 2 ,^) = F^(Q 2 , -a;). (2.27) 

In pQCD calculations higher-order corrections to various physical processes, are, in gen- 
eral, large and in order to improve the convergence of the corresponding perturbative series, 
different methods have been proposed. In exclusive processes - considering the pion elec- 
tromagnetic form factor as being a prominent example - the next-to- leading order 0(a s ) 
correction contains logarithmic terms ~ \n(Q 2 xy / /z 2 ^) (or ~ \ia(Q 2 xy/ 0, IS], which 
can be entirely or partly eliminated by a proper choice of the renormalization scale \j? R (]!%)■ 
This can be achieved by taking as a renormalization scale \j? R = Q 2 xy [p 2 R = Q 2 xy] or 
fi 2 R = Q 2 x/2 {p 2 R = Q 2 x/2}. The renormalization scale enters into the pQCD expression not 
only via logarithmic terms, but also through the argument of the running coupling a s (fi 2 R ). 
In order to calculate the amplitude of an exclusive process, an integration in corresponding 
integrals over the longitudinal momentum fractions of the quarks and gluons in the involved 
hadrons has to be carried out. Thus, in the pion electromagnetic FF calculations, the in- 
tegration over the variables x and y has to be performed. Traditionally, to avoid problems 
associated with singularities of a s (fi 2 R ), one "freezes" the running coupling by replacing x 
and y, e.g., by their mean values x — >< x >= 1/2, y — » < y >= 1/2 in the argument 
of a s {^ 2 R ) and performs then the calculation with a s (Q 2 /A) |44j. Recently, the RC method 
for the calculation of various exclusive processes with a s (/•*#) was proposed. Within this 
framework, one expresses the running coupling a s (Q 2 x) in terms of a s (Q 2 ) by employing 
the renormalization-group equation and performs the integration over x (y) using the prin- 
cipal value prescription. This allows one to estimate power-suppressed corrections to the 
process under consideration towards the endpoint region of phase space, thus improving the 
agreement of QCD theoretical predictions with experimental data. An alternative option 
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to calculate power-behaved contributions with their coefficients would be to use the ana- 

Next-to-leading order 
In general, the renor- 



lytic perturbation theory |33j along the lines proposed in 



corrections are known for only a few exclusive processes |44l . 1451 14 
malization scale n 2 R may be taken equal or proportional to the square of the four-momentum 
q 2 of the virtual parton(s) mediating the strong interaction in corresponding leading-order 
Feynman diagrams. 

For the //-meson on-shell gluon transition (i.e., for u = ±1), we adopt in this paper the 
choice 



/4 



Q 2 x, 



/4 



Q 2 x 



(2.28) 



The renormalization scales ()2.28|) are equal to the momentum squared \q 2 \ of the virtual 
partons (gluon or quark) in the corresponding Feynman diagrams. In the general case 
(lo 7^ ±1), the absolute value of the square of the four-momenta q 2 of the virtual partons 
depends on both the total gluon virtuality Q 2 and the asymmetry parameter u. However, 
for the sake of simplifying the calculations and to avoid problems related to the appearance 
of two parameters Q 2 and u in the argument of a 8 , we shall use a renormalization scale of 
the form given in Eq. ()2.28|) . This choice is justified from the physical point of view because 
those parts of the scales ~ Q 2 x, ~ Q 2 x are exactly responsible for the power corrections 



i/Q 



2 11 



n 



1, 2, .. to the form factor Q F v i g * g *(Q ,u) which we are going to compute. 



The next problem then is how to generalize Eqs. (j2.19|) and ()2.20|) in such a way so that 
the corresponding hard-scattering amplitudes T^) and the quark and gluon components of 

For this purpose, we symmetrize the functions 
~p 2 R to obtain 



the vertex function will obey Eqs. 
T*g } (a;,Q 2 ,u;,/4(7Z|)) by exchanging /j? R 



T?(x,Q 2 ,u) = - T?(x,Q 2 ,u,ti 2 R )+T?(x,Q 2 ,u;,iI 2 R 



2tt a s (Q 2 x) + a s (Q 2 x) 
3Q 2 (1 + u)x + (1 - u)x' 

T$(x, Q 2 , u, jf R ) + T$(x, Q 2 , u, n 2 R ) 

2tt a s (Q 2 x) +a s (Q 2 x) 
3Q 2 (l + u)x+ (1 -u)x 

T? (x, Q 2 , u, //|) + Tf (x, Q 2 , u, Ji 2 R ) 
(1 + uj)x + (1 — to)x 



and 



T«(x,Q 2 ,u) 



Tf(x,Q 2 ,u) 



71 



71 



2Q 2 n } 



2Q 2 n f 
T$(x,Q 2 ,uj 



a s (Q 2 x) + a s (Q 2 x 
1 



LU [(1 + u)x + (1 — Oj)x] 

Ti{x,Q 2 ) u ) jl 2 R )+Ti{x 1 Q 2 ^^ 2 R ) 

2 _ A ] (1 + Oj)x + (1 — U))x 

uj [(1 + u)x + (1 — uj)x] 



a s (Q x) + a s (Q x 



[2.29) 



(2.30) 
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In the "frozen" coupling-constant approximation Q 2 x, Q 2 x — > Q 2 , Eqs. (|2.29|) and ()2.30|) 
coincide with expressions ()2.19|) and ()2.20|) . It is also not difficult to verify that the hard- 
scattering amplitudes T?$(x, Q 2 , u) (l2~2H|) . (I2~3U1) satisfy Eqs. (EHHj) and RFM . As a result, 
Eqs. fl2.25|) - ()2.27|) remain valid also within the RC method. 



C. Quark and gluon components of the 7/ meson distribution amplitude 

The important input information needed for the computation of the form factor 
F v 'g* g *(Q 2 , o>) are the DA's of the quark and gluon components of the rj' meson, namely, 
the functions <fr q (x,Q 2 ) and <ft 9 (x,Q 2 ). In general, a meson DA is a function containing all 
nonperturbative, long-distance effects, which cannot be calculated by employing perturba- 
tive QCD tools. Its dependence on x (or its shape) has to be deduced from experimental 
data or found using some nonperturbative methods, for example, QCD sum-rules (see, e.g., 
|47l|). In contrast, the evolution of the DA with the factorization scale Q 2 is governed by 
pQCD. 

The evolution equation for the DA of a flavor singlet pseudoscalar meson was derived 

and solved in Rets. BE3II3- 

It turned out that due to mixing of the quark- ant iquark and 



two-gluon components of the meson DA, the evolution equation has a (2 x 2) matrix form. 
The anomalous dimensions matrix, which enters into the evolution equation at the one-loop 
order, has the following components (UHH EH (see also Ref. [iil ]) 



7, 



m 



c h 



(n + l)(n + 2) 



n+l i 
3=1 J 



' ' 99 iVc 



Yin 



la 



f 



7 



c f 



Po 8 

N c (n+l)(n + 2) 

n(n + 3) 



n+l i 
3=1 J 



[2.31] 



,q9 (n + l)(n + 2)' lgq r 3(n + l)(n + 2)' 

where N c = 3 and Cp = 4/3 is the group theoretical factor for SU C (3). The anomalous 
dimensions matrix has the eigenvalues 

„ 1 



7± = 2 



199 + ll ± y/(%, ~ 199) 2 + 



(2.32) 



The solutions of the evolution equation for the quark and gluon components of the DA has, 
in general, the form 



(x, Q 2 ) = 6Cxx 1 1 + jr 



Bl 



k=2A.. 



00 



+ p 9 B 9 

1 rnn 



00 



CT(x-x) 



and 



(2.33) 



[x, Q 2 ) = Cxx ^2 



PlB 



k=2A.. 



MQ 2 ). 



00 



+ B 9 

1 n 



MQ 2 ), 



00 



cj_i(x 



X) 



(2.34) 
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with the constant C being defined as 

C = V2f q sin P + f s cos < 



In Eqs. (jUHHl) and (ESlC'n /2 (^) and C£/ 2 (z) are Gegenbauer polynomials, calculable using 
the recurrence formula [49] 

(n + l)C» +1 (z) = 2(n + v)zCl{z) - (n + 2u - lJC^z), (2.35) 

C£(z) = 1, C[{z) = 2vz. 

The parameters and are determined in terms of the anomalous dimensions matrix 
elements 

p i=q1± r M p 9 t (2.36) 

'99 '99 

In Eqs. (I2.33j) and (J2.34)) /io = 1 GeV is the normalization point, at which the free input 
parameters B%, B 9 in the DA's 4> q (x,Q 2 ) and 4> 9 (x,Q 2 ) have to be fixed. Exactly these 
parameters determine the shape of the DA's. 

For all n > 2, the eigenvalues 7™ < and their absolute values increase with n. Hence, 
in the asymptotic limit Q 2 — > 00, one has 



00 



MQ 2 ) 



ln(g 2 /A 2 ) a °~ ^0 



and therefore only the quark component of the 7/ meson DA survives, evolving to its asymp- 
totic form, whereas the gluon DA in this limit vanishes, 

<? (x,Q 2 ) Q ^° QCxx, <p 9 {x,Q 2 ) Q2 -^ 0. (2.37) 

In this work, we shall use the r^'-meson DA that contains only the first non-asymptotic 
terms. In other words, we suggest that in Eqs. (j2.33|) and ()2.34|) B\ 7^ 0, £>f 7^ and 
B q = B 9 = for all n > 4. The numerical values of the relevant parameters for rif = 3 are 

2 _ 2 _•! I 2 ^ 2 _ o 

7gg — g ) Igg — LL i Igq ~ ^99 ~ ' 

2 48 2 101 - 16 1 

7 2 ~ , 7 2 ~ , pi ~ — , pi ~ . (2.38) 

Taking into account Eq. (|2.38j) and the expressions for the required Gegenbauer polynomials 



C\ [x — x) — - 5(x — x) 2 — 1 = 6 (1 — 5xx) , C^ /Z (x — x) — 5(x — x), 

we finally recast the rj meson quark and gluon DA's into the more suitable (for our purposes) 
forms 

(f)i(x, Q 2 ) = 6Cxx \l + A(Q 2 ) - 5A(Q 2 )xx] , <j) 9 (x, Q 2 ) = Cxx(x - x)B(Q 2 ). (2.39) 



5/2, 
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In these expressions the functions A(Q 2 ) and B(Q 2 ) are defined by 

A{Q 2 ) - 




B(Q 2 ) = 16BU^y-\ +5j Bfp^ . (2.40) 



The rj' meson quark and gluon DA's for rif = 4 are given in the Appendix. Equations (|2.39|) . 
(I2.4U)) . and (|A.1|) contain all necessary information about the singlet part of the 7/ meson 
DA. 



III. RUNNING COUPLING METHOD AND POWER SUPPRESSED CORREC- 
TIONS 



To compute the //-meson virtual-gluon transition form factor F v i g * g * (Q 2 , u) , we have to 
perform in Eqs. (j2.14|) and (j2.15|) the integration over x by inserting the explicit expressions 
for the hard-scattering amplitudes T^\x, Q 2 , u) and the 7/-meson quark and gluon DA's, 
and retain the x-dependence of the coupling a s (Q 2 x) [a s (Q 2 x)}. Such calculations lead to 
divergent integrals because the running coupling a s (Q 2 x) [a s (Q 2 x)} suffers from an infrared 
singularity in the limit x — > [x — > 1]. This means that in order to carry out computations 
with the running coupling, some procedure for its regularization in the end-point x — > 0, 1 
regions has to be adopted. 

As a first step in this direction, we express the running coupling a s (Q 2 x) in terms of 
a s (Q 2 ), employing the renormalizat ion-group equation 



da s (Q 2 x) 
d In x 



A 

47T 



a s (Q 2 x) 



8tt 2 



a s (Q 2 x) 



(3.1) 



The solution of Eq. (|3.1|) . obtained by keeping the leading (a s \nx) k and next-to-leading 



a s (a s lnx) fe 1 powers of Ins, reads [5 



a s (Q 2 x) 



« S (Q 2 



l + \nx/t 



1 - 



a s (g 2 )Aln[l + lnx/t] 
2vr/3 l + lnx/t 



(3.2) 



where a s {Q 2 ) is the one-loop QCD coupling constant and t = 4n/ (3 a s (Q 2 ). 

Inserting Eq. ()3.2j) into the expressions for the hard-scattering amplitudes and subse- 
quently T12 \x, Q 2 , u) into Eqs. ()2.14|) and (|2.15|) . we obtain integrals which are still di- 
vergent, but which now can be calculated using existing methods. One of these methods, 
applied in Ref. [3cj for the calculation of the electromagnetic pion form factor, allows us 
to find the quantity under consideration as a perturbative series in a s (Q 2 ) with factorially 
growing coefficients C n ~ (n — 1)!. Similar series with coefficients C n ~ (n — 1)! may be 
found also for the form factor Q 2 F r] i g * g * (Q 2 , u), i.e., 



Q Frf g*g* (Q 



00 



E 

n=l 



a s {Q 2 ) 
An 



(3.3) 



But a perturbative QCD series with factorially growing coefficients is a signal for the IR 
renormalon nature of the divergences in Eq. (j3.3J) . The convergence radius of the series ()3.3|) 
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is zero and its resummation should be performed by employing the Borel integral technique. 
First, one has to find the Borel transform B[Q 2 F r) i g * g *\{u) of the corresponding series [5 1| 



B 



Q 2 F n 



9 9 



oo u n-l 



E 

71=1 



n 



1 ! 



:C n . 



(3.4) 



Because the coefficients of the series f!3.3|) behave like C n ~ (n— 1)!, the Borel transform (|3.4|) 
contains poles located at the positive u axis of the Borel plane. In other words, the divergence 
of the series f!3.3|) owing to Eq. ()3.4|) has been transformed into pole-like singularities of the 
function B[Q 2 F v i g * g *](u) ■ These poles are the footprints of the IR renormalons. 

Now in order to define the sum ()3.3|) or to find the resummed expression for the form 
factor, one has to invert B[Q 2 F v > g * g *](u) to get 



[Q 2 F v , g * g *(Q 2 ,u) 



P.V. 



du exp 



47T 



(3 a s (Q 2 



-u 



B 



[U 



(3.5) 



and remove the IR renormalon divergences by the principal value prescription. 

The procedure described above, is straightforward but, at the same time, tedious. For- 
tunately, these intermediate operations can be bypassed by introducing the inverse Laplace 
transforms of the functions in Eq. ()3.2j) [Uj], viz., 



1 



and 



(t+zy 

\n(t + z) 
(t + z) 2 



r 



i f°° 

— / du exp[— u(t + z)\u v ~ , 
\v) Jo 



Rev > 



/■OO 

/ duexpl— u(t + z)](l — 7b — In -u)ii, 
Jo 



(3.6) 



(3.7) 



where T(z) is the Gamma function, — 0.577216 is the Euler constant, and z = lnx [or 
z — Inx]. 

Then using Eqs. (Q, (jSSl) and (jSZ3) for a s (Q 2 x), we get @ 



ot s (Q 2 x) = a s {Q 2 )t / duexp(—ut)x~ u R(u,t) = — / duexp(—ut)x~ u R(u,t) (3.8) 
Jo Pq Jo 

with the function R(u, t) defined as 

R(u,t) = 1 - t^u(1 -1e -lnt-lnu). (3.9) 

Po 

Employing Eq. ()3.8|) for a s (Q 2 x) and carrying out the integrations over x, we obtain (see, 
the next Section) the quark and gluon components of the form factor F^^g* (Q 2 , uj) directly 
in the Borel resummed form ()3.5j) . 

The resummed vertex function [Q 2 F v i g * g * (Q 2 , uj)] res - except for the //-meson on-shell- 
gluon transition and the uj = case - contains an infinite number of IR renormalon poles 
located at Uq = n, where n is a positive integer that, in general, is given as a sum of a series, 
i.e., 



CO 



Qtu n2 ^ \- 4?r D TA r e- ut R(u,t)du 

S(k ,Q ,uj) ~ 2^ ir P - V - / C n (u). (3.10) 

n =k Po Jo n-u 
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One appreciates that all information on the process asymmetry parameter u is accumulated 
in the coefficient functions C n {uo). 

To avoid inessential technical details and to make the presentation as transparent as 
possible, let us for the time being neglect the non-leading term ~ a 2 in Eq. (|3.2j) and make 
the replacement R(u,t) — > 1 in Eq. (|3.1U|) . Then, the integrals in Eq. (|3.1U|) take, after 
simple manipulations, the form 



f RV - 



e- ut du Air li(X r 



X = Q 2 /A 2 



Jo n — u (5q \ n 
with the logarithmic integral defined by 



(3.11) 



The function is expressible in terms of the function Ei(z) |4 

/z(0 = -^i(-lnO 

where 



(3.12) 



(3.13) 



Ei{z) = L 



e' l r l dt 



The latter expression can, for \z\ — > oo, be expanded over inverse powers of z 

M 



EAz) 



z-\- z 



E 



ml 



Using Eqs. (l3~T3]) and flTHJ), we find 

t M 



x - mi 



o ("^) m 



+ o(l 



-(M+l) 



M 



(3.14) 



It is not difficult to see that for Q 2 ^> A 2 



4n li(X n ) 
To A« 



M 



4-E 

m=l 



a s (Q 2 

47T 



1 ^ ml 
"^R^ (nlnO ? 



(3.15) 



Po 



,-i(m-l)\ 



71" 



(3.16) 



Choosing in Eq. ()3.16|) M — > oo and comparing the obtained result with Eq. ()3.3|) . we 
conclude that each term in the sum ()3.10|) at fixed n is the resummed Borel expression of 
the series (|3.16|) . 

But we can look at the integrals in Eq. ()3.1()j) from another perspective. Namely, it is 
easy to show that 



47r r°° e ut du 
Po Jo n — u 
where f2n(Q) are the moment integrals 

f P (Q) 



i i 

a s (Q 2 x)x n ~ 1 dx = -f 2n (Q), 
n 



Qp Jo 



Q 



dkk p ~ l a s (k 2 ). 



(3.17) 



(3.18) 
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The integrals f p (Q) can be calculated using the IR matching scheme [5 

p N 



fp{Q)=(%) / P M + «s(Q 2 )E 



n=0 L 



A) 

a s {Q 



2np 



[n!-r(n + l, pHQ/fu))}, (3.19) 



where fij is the infrared matching scale and T(n + 1, z) is the incomplete Gamma function 



OC 



r (n + 1, z) = / e _t t n dt. (3.20) 

In Eq. (I3.19|) f p (fii) are phenomenological parameters, which represent the weighted 
average of a s (k 2 ) over the IR region < k < fii and act at the same time as infrared 
regulators of the r.h.s. of Eq. ()3.17j) . The first term on the r.h.s. of Eq. ()3.19j) is a power- 
suppressed contribution to f p (Q), which cannot be calculated within pQCD, whereas the 
second term is the perturbatively calculable part of the function f p (Q). Stated differently, the 
infrared matching scheme allows one to estimate power corrections to the moment integrals 
by explicitly dissecting them from the full expression and introducing new nonperturbative 
parameters The values of these parameters have to be extracted form experimental 

data. 

In order to check the accuracy of both the RC method (IR renormalon calculus) and the 
infrared matching scheme, a Landau-pole free model for a s (Q 2 ) was introduced in Ref. 
yielding 



« S (Q 2 ) = ^ 



1 +4A 

+ 125- 



(3.21) 



In A (1-A)(4 + A) 4 _ 

The model expression (|3.21j) was used for the exact calculation of the moment integrals 
/p Xact (Q). The IR renormalon and the IR matching scheme results were compared with 
/p Xact (Q) concluding that for N > 2 the IR matching scheme leads to a good agreement 
with the exact values of f p (Q). As regards the calculation of the integrals f p (Q), when 
employing the RC method with the principal value prescription, a good agreement with 
/| xact (Q) was found for p = 2, while for p = 1 a small deviation from /f xact (Q) in the region 
of Q 2 ~ a few GeV 2 was observed (the lowest moment integral entering into our expressions 
is fziQ))- But even for p = 1, the agreement with /f xact (<5) within the error bars is quite 
satisfactory in view of the uncertainties produced by the principal value prescription itself. 
Therefore, we can conclude that the RC method and the IR matching scheme lead for the 
moment integrals f p (Q) to approximately similar numerical results. It is remarkable that in 
both cases an enhancement of the IR renormalon and matching-scheme results relative to 
the leading-order perturbative prediction was found. 

After presenting our analysis, it becomes evident that the Borel resummation technique 
enables us to estimate with a rather good accuracy power-behaved corrections to the integrals 
in Eq. (|3.1U|) and hence to the vertex function F v > g * g *(Q 2 , to). Of course, one can employ 
the IR matching scheme and Eq. (J3.19j) to determine power corrections to the form factor 
in explicit form. But as we have noted above, in general, the form factor Q 2 F rj i g * g * (Q 2 , to) 
contains an infinite number of IR renormalon poles or stated equivalently, an infinite number 
of f2n(Q) is required to calculate sums like the one presented in Eq. ()3.10|) . In real numerical 
computations, in order to reach a good approximation of the series (J3.10)) . at least a sum 
of the first 15 4- 20 terms is needed, meaning the involvement of /Wi (/-£/) [n = 15 4 20) 
nonperturbative parameters. But from the experimental data, only the first few moments 
fp(f^i) could be extracted [H^, E^, 0, Values for / x (2 GeV) (in the literature the 
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notation /i(/ir) = ao(/ij) is widely used) range from ao(2 GeV) = 0.435 ± 0.021 [54j and 
« (2 GeV) = 0.513t° fi Q to a (2 GeV) = 0.597±8;g?g @. Here, we write down only 
sample results, obtained in these works, in order to demonstrate the experimental situation 
with a (2 GeV) = f x (2 GeV). For / 2 (2 GeV) the value / 2 (2 GeV) ~ 0.5 was deduced 



53, 581 from the data on structure functions 57 



In the framework of the RC method, we estimate the same power corrections to a physical 
quantity, but here we do not need additional information on / p (///). Moreover, this method 
gives us the possibility to calculate values of the nonperturbative parameters f p (ni) by 
means of the formula 

fM = ^TT-rA A = /i 2 /A 2 (3.22) 

that can be easily found by comparing Eqs. (J3.11|) and (J3.17|) . The calculated values of the 
first few even moments (rif = 3, A = 0.3 GeV) are 

/ 2 (2 GeV) ~ 0.535, f 4 (2 GeV) ~ 0.45, f 6 (2 GeV) ~ 0.41. (3.23) 

To compare our predictions with those computed via the model a s (Q 2 ) (Eq. (|3.21j0 . we 
choose rif = 3 and A = 0.25 GeV. The corresponding values are shown below 



/ 2 RC (2 GeV) ~ 0.479, f 4 RC (2 GeV) ~ 0.393, ., 
/ 2 mod (2 GeV) ~ 0.450, ff od (2 GeV) ~ 0.388 °~ ' 



One observes that the parameters f p {fii) of the RC method are in agreement with both the 
experimental results and the model calculations. 

We have noted above that the principal value prescription, adopted in this work in order 
to regularize divergent integrals (see, Eqs. (|3.5)L (|3.10Jl ) generates power-suppressed (higher- 
twist) uncertainties 

~L, N *—Q3T> ( 3 - 25 ) 

where {§ q (Q 2 )} are calculable functions entirely fixed by the residues of the Borel transform 
B[Q 2 F Mg *g*]{u) at the poles q = u and {N q } are arbitrary constants to be fixed from 
experimental data. In the case of the 7/7 electromagnetic transition form factor Q 2 F r] r 1 (Q 2 ), 
the uncertainties (|3.25J) were estimated in Ref . Q . For the 7/ meson asymptotic DA and for 
the parameters N± = N2 = 1 and N\ = V 2 = —1, it was found that they do not exceed ±15% 
of the 7/7 form factor. Because the RC method allows one to evaluate power corrections 
to a physical quantity, in general, and for the r/x transition form factor it has already 
provided a good agreement with the CLEO data |6(, one has to introduce the constraint 
1-^1,2 1 ^ 1 in order to retain this agreement. In the present work we shall estimate in our 
numerical analysis the higher-twist uncertainties produced by the principal value prescription 
by choosing the constants {N q } — ±1. It will be demonstrated that the ±15% bounds are 
valid also for the rj'g, rfg* transition FF's. Inclusion of higher-twist corrections, arising 
from the two particle higher-twist and higher Fock states DA's of the rj' meson, may, in 
principle, further improve this scheme. These refinements are, however, beyond the scope 
of the present investigation. 
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IV. THE FORM FACTOR F v , g * s * (Q 2 , u) WITHIN THE RC METHOD 

In this section we calculate the quark and gluon components of the vertex function 
F v > g * g *(Q 2 , a>) within the RC method. We also present our results for F v / gg *(Q 2 ,uj = ±1), 
Fjf g * g * (Q 2 , oo = 0) and the asymptotic limit for the form factor. 



A. Quark component F^, g „ g „ (Q 2 , lo) of the vertex function 

To calculate the quark component of the vertex function Fl g * g „ (Q 2 , u), we use Eq. (I2.25j) 
for F^, g * g *(Q 2 , uj), Eq. (|2.29|) for the hard scattering amplitude T±(x, Q 2 , u), and the expres- 
sion for a s (Q 2 x) given in Eq. (|3.8jl . Employing the quark component of the 7/-meson DA, 
(f) q (x, Q 2 ), we get 4 

Po J o 



X 



x u xdx 



+ 



o (1 + uj)x + (1 — oj)x Jo (1 + u)x + (1 — lo)x 



xx u dx 



160n 2 CA(Q 2 

Po 



POO 

■ / due~ ut R{u,t) 
Jo 



(1 + to)x + (1 — Uj)x 



+ 



1 X 2 X 2 u^ x 



10 (1 + u)x + (1 — U))X 

Using the fact that the integral jH^] 

f 1 x a - l x^-\l - xrVdx = B(a, p) 2 F 1 (-7, a; a + (3; r) 
Jo 

is expressible in terms of the Gauss hypergeometric function 

00 (a) t (b) t z» 



2 F l(a ,6;c;z)=g ^ fc . 



we obtain 



Q 2 F q vl g,g,{Q\0O) 



32n 2 C[l + A(Q 2 )} 



/ due- ut R{u,t)B(2-u,2) 
Jo 



2u \ ( 2cu 

2 F 1 {l,2;A-u;^—)+ 2 F 1 l,2- U ;4-u; 



1 +UJ 



160tt 2 CA(Q 2 ) r°° 



/•OO 

/ due~ ut R(u,t)B(3-u,3) 
Jo 



1+LU 



Fx 1,3:6 -u) 



2oo 

1+UJ 



(4.1) 



(4.2) 



(4.3) 



4 In what follows integrals over the variable u have to be understood in the sense of the Cauchy principal 
value. 
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+2-Fi 1 1, 3 - u: 6 — u: 



(4-4) 

In Eqs. (14 .2)1 and (J4.3|) . 5(x, y) and (a) n are the Beta function and Pochhammer symbols, 
respectively, defined in terms of the Gamma function T(z) 



r(x+y)' 1 "' I» ' 

From Eq. (|4.4|) a simple expression for the r/-meson on-shell-gluon transition form factor 
F^, gg *(Q 2 , uj = ±1) can be found 

Q 2 F q n , gg *{Q 2 ^ = ±1) = - 167r2C[1 ^ A(g2)] J™ due-^R(u,t) [5(1,2 - u) 

+5(2, 1 - u)] + 807r2< "l A(g2) /°° due- ut R(u, t) [5(3, 2 - u ) + 5(2, 3 - u)} , (4.5) 

Po Jo 

where we have used the equality 

^/ , \ r(c)r(c — a — b) 

2 F 1 (a,b;c;l) = J 1 -f 4.6 

r(c — a)T(c — b) 

It is worth noting that Eq. (|4.5|) can be obtained from Eq. (J2.25j) by employing the uo — > ±1 
limits of the hard scattering amplitudes (|2.29|) 



T?{x, Q\ u = ±l)+ T$(x, Q 2 , u = ±1) 



7T 



3Q 2 



a e (Q x) + a s (<5 x 



1 1 
- + = 

X X 



(4.7) 



In the case of gluons with equal virtualities Q\ = Q\ (and hence in the uo = case), the 
form factor can be found employing the expression for the hard-scattering amplitude 



2?(x, Q 2 , w = 0) + T|(x, Q 2 , = 0) 
Calculation leads to the simple expression 



An 



3Q 2 



a s (Q 2 x) + a s (Q 2 x) 



(4.8) 



Q*F} g . g .{Q\u> = 0) 



64n 2 C(l + A(Q 2 )) 
Po 



/ due- ut R{u,t)B(2-u,2) 
Jo 



+ 



320tt 2 C7L(Q 2 ) /-°° 



0o 



poo 

/ ciue- u *5(M,t)5(3-u,3). 



(4.9) 



Let us note that Eq. (|4.9jl can be deduced from Eq. (|4.4|) in the limit a; — > by taking into 
account that [HO ] 

a F 1 (a,6;c;0) = l. 

The next important problem to be solved within the framework of the RC method, is 
to reveal the IR renormalon poles in Eqs. ([4.4)1 . (|4.5jl and ()4.9|) because without such a 
clarification all these expressions would have merely a formal character. 

We start from the simple case, i.e., from Eq. ()4.5jl . and get 



5(1, 2 -m) + 5(2,1 -n) 



r(i)r(2-«) r(2)r(i-«) 



r(3 - u) 



r(3-«) 
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1 

+ 



2-u (l-u)(2-u) l-u 
and 

r(3)r(2-«) r(2)r(3-«) 



5(3, 2 - u) + 5(2, 3 - u) 



+ 



r(5-«) r(5-«) 
1 1 i 



(2-w)(3-w)(4-«) (3-«)(4-«) 2-u 3-u" 
In deriving these expressions we used the following property of the Gamma function 

zF(z) = T(z + l). 

Hence we have 

' yy po «/o 1 — u 

+ ? 2! f^y- ,1 1^ 

p jo \2 — u o — u/ 

As one sees, the structure of the IR renormalon poles is very simple; we encounter only a 
finite number of single poles located at Uq = 1, 2 and 3. 
In the same manner we get from Eq. ([4.9]) 



/5o Jo v ' y (2-w)(3-w) 

+ 640,W) r 1 

(3 Jo (3 - u)(4 - u)(5 - it) 

In this case the single IR renormalon poles are located at the points Uq = 2, 3, 4 and 5. 

In order to get the IR renormalon structure of the integrands in Eq. (|4.4|) . we have to 
expand the hypergeometric function 2 Fi(a,b; c; z) in powers of z (see Eq. ([4.3]) ). where the 
condition \z\ < 1 must hold. But the argument of the function 2 Fi(a,b;c;2uj/(l + uj)) 
satisfies this requirement only in the region uj G (0,1). In the region u G (—1,0) an 
expression obtained from (J4.4j) by means of a simple transformation (see, (J4.17)) below) 
has to be used because in this case the argument of the hypergeometric function becomes 
equal to 2uj/(uj — 1) < 1, obeying in the region uj G (—1,0) the required constraint. But 
regardless of the expansion region, we obtain in both cases the same IR renormalon structure. 
Adding to this argumentation the evident fact that the vertex function is symmetric under 
the replacement uj <-> —uj (cf. Eq. ([2.27]) ). we can restrict the study of the form factor 
Q 2 F5 g * g *(Q 2 , uj) to the region uj G (0, 1). Then, we can expand the hypergeometric functions 
2-F\(ct, b; c; 2u/(l + uj)) in the region uj G (0, 1), via Eq. ([4.3]) . For example, for one of these 
functions, we get 

B(2<2 , ( U/1H ?j .^_ r(2)r(2-n) - (l) fc (2) fc /3 fc 
3(2, 2 - u)^ (1, 2, 4 - u, 0) - r(4 _ M) £ (J3^X! 

" r(fe + 2)r(2-^) ~ 

"to r (fc + 4-n) * "fro ( + ' ^ 
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with 

a 2w 



1 + UJ 

The remaining terms in Eq. ()4.4)1 can be treated in the same manner and as a result we 
obtain 



oo 



X 



[B{2 - u, k + 2) + 5(2, k + 2-u)]/3 



k=0 



+ ^ l ; ' I" due- ut R{u, t) Y [5(3 - u, k + 3) + 5(3, k + 3 — u)] f3 k . (4.12) 

The IR renormalon structure of the integrands in Eq. ()4.12|) is quite clear now. In fact, we 
can write the Beta functions entering Eq. (j4.12|) in the following form 



, , . r(2)r(fc + 2-«) 

B(2,k + 2-u) - 



T(k + A-u) (k + 2-u)(k + 3-u) 
and correspondingly 

5(2 — u,k + 2) — w T(y \ +2 }i s, 

y ' ' (2-u)(3-u)...(k + 3-uY 

B(3,k + 3-u) = 
B(3-u,k + 3) 



(k + 3 - u)(k + 4 - u)(k + 5 - u) 
T(k + 3) 



(3-u)(A-u)...(k + 5-u)' 

Here we have an infinite number of single IR renormalon poles located at the points Uq = 
k + 2, k + 3; uq — 2, 3, . . . k + 3; u = k + 3, k + 4, k + 5 and Uq — 3, 4, . . . k + 5, respectively. 

The last question to be answered is, whether one can use our results, obtained in the 
context of the RC method, in the limit Q 2 —>■ oo in order to regain the asymptotic form of 
the form factor F^, gtgt (Q 2 — > oo,u). It is clear that regardless of the methods employed and 
the approximations done, in the limit Q 2 — ► oo the form factor F^, g * gt (Q 2 ,ui) must reach 
its asymptotic form. This is true, of course, for our computations, as we estimate power- 
suppressed corrections to the form factor Q 2 F^, g , g , (Q 2 , u>), which become important in a 
region of Q 2 ~ of a few GeV 2 , but vanish in the asymptotic limit. As we have emphasized 
in Sect. II, in the asymptotic limit the gluon DA of the 7/-meson satisfies <ft 9 (x, Q 2 ) — > and 
hence 

The DA of the quark component <f) q (x,Q 2 ) of the r( meson evolves for Q 2 —>■ oo to the 
asymptotic DA ()2.37j) and all non-asymptotic terms in <fi q {x,Q 2 ) proportional to C^J 2 {x — 
x), n > (in our case ~ A(Q 2 )) vanish. Therefore, the results, which we shall obtain here, 
describe not only the asymptotic limit of Q 2 F q lgt , gt ,(Q 2 ,uj), but also the asymptotic limit of 
the vertex function Q 2 F v i g * g *(Q 2 ,u) itself. 

In the limit Q 2 — > oo the asymmetry parameter can take values to — > ±1 (if we pass to 
the limit Q 2 — >• oo at fixed Q\ or Q\), uj = (Q\ = Q\ and Q 2 — > oo) or uj ^ ±1, (if we 
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take the limit Q 2 — > oo at fixed u;). We consider here all possibilities: a) Q 2 — > oo, c<j — > ±1, 
b) Q 2 — > oo, u; = and c) Q 2 — > oo, ±1, 0. In the limit Q 2 — > oo, we also take into 
account that the second term in the expansion a s (Q 2 x) ()3.8|) has to be neglected. In other 
words, in the limit Q 2 — > oo, we find 



/■OO /"OO 

/ e~ ut R{u,t)du^ / e" 

JO JO 



Iff 



(4.13) 



We begin from the simpler case a). In Eq. (j4.1U|) we have already obtained the desired 
limit, but A(Q 2 ) ^ 0. Taking into account Eq. (l4~T3l and A(Q 2 ) -> 0, we get 



Q 2 F v , gg *(Q 2 ,LU = ±l) 



lQn 2 C e~ ut du 16tt 2 C li(X) 



Po Jo 1 - u 



00 



A 



It is easy to show that using only the leading term in the expansion of li(X)/X (see, Eq. 
(J3HJ)), we find 

Q 2 F ri , gg *(Q 2 ,u = ±1) Q -^° -4nCa s (Q 2 ). (4.14) 
The limit Q 2 —>■ oo, uj = can be analyzed by similar means. Thus, from Eq. (|4.11|) we get 



Q 2 F v , g * g *(Q 2 ,u = 0) 



64vr 2 C 



e du 



64n 2 C 
Jo 



Po Jo (2-u)(3-u) 

u(x 2 ) ii(x s y 



x 2 



A 3 



which in the limit under consideration simplifies to 

ga-oo 87rCa s (Q 2 ) 



g 2 F r?w (g 2 ,w = o) 



(4.15) 



Now let us consider the more interesting case c). Then, from Eqs. (J4.4j) and (J4.13|) . we 
obtain 2 

Q 2 F v ' g * g *(Q 2 ,uj) Q — 



A)(l + w) io 



due~ ut B 2, 2 - u) 



x 



Fi l,2;4-u; 



2cj 
1+w 



+ 2 Fi l,2-u;4-u; 



2cj 
1 + 



(4.16) 



As an example, we analyze the second term in Eq. (|4.1fij) (see also Eq. (|4.12jl ) 

POO POO °° 

/ due~ ut B(2,2 -u) 2 F 1 (1,2 -u; 4 - u; f3) = due~ ut V B(2, k + 2 - u)(3 
Jo Jo f-i 



x fOO 

Y / T(2)p k due 
k=o Jo 



Lit 



1 



1 



k + 2 — u k + 3 — u 



£ r (2) 

fc=0 



u(x k+2 ) h(x k+3 ) 



x k + 2 



In the considered limit, one finds 

^(A fc + 2 ) U(X k+3 ) 1 



X k+2 



X k+3 



1 



1 



lnAVA; + 2 k + 3 



In A (fc + 2)(fc + 3)' 
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Then, the result reads 

~h(\ k+2 ) u(\ k + 3 ) 



£r(2) 

k=0 



\k+2 



3 k -> — V 



£ B(2, k + 2)(3 k = t-^B{2, 2)2^(1, 2; 4; /?). 



lnA^ ~ In A 



The same method can be applied to the first function in Eq. (|4.16|) . But before doing that, 
it is instructive to employ the transformation 



2 



Ftia, b, c; z) = (l- z)-\F x [a, c-b,c; — ^) . (4.17) 



After performing all these operations we find 



Q 2 F, W (Q 2 , u) - ™* 2 C -L B (2, 2)^(1, 2; 4; (3) 

" ' 2^(1, 2; 4; /?). (4.18) 



3(l+a;) 

Taking into account the expression for the hypergeometric function 2-^1(1, 2; 4; /?) in terms 
of elementary ones 



2^(1, 2; 4; (3) = ^ [(3(2 -(3) + 2(1 - (3) ln(l - 0)] , (4.19) 



we finally obtain 



Q 2 F, lg , g ,{Q 2 ,u) -4vrCa s (Q 2 )/M, /(c) = -I ( 1 + hi i— - ) . (4.20) 

In our argumentations we have tacitly assumed that uj G (0, 1). But Eq. (j4.20j) holds for all 
values of uj ^ ±1, 0, which is evident from the equality f{uj) = f(—u). 

Equations (|4.14|) . ()4.15|) and ()4.20|) can be obtained in the standard HSA using the 
corresponding hard-scattering amplitudes ([4.7)1 . (|4.8|) and (|2.19|) . The analysis carried out 
so far proves the correctness of the RC method leading to the expected expressions for the 
form factor Q 2 F v > g * g * (Q 2 , uj) at Q 2 — > oo and demonstrating at the same time the consistency 
of the symmetrization of the hard-scattering amplitudes ()2.29|) . ()2.30|) . 



B. Contribution of the gluon component of the r/-meson to the form factor 

F^/g.g* (Q 2 , u) 

The contribution of the gluon component of the ?/-meson to the form factor F v > g * g * (Q 2 , uj) 
can be computed using the methods described in the previous subsection. 

Using Eq. ()2.30|) for the hard-scattering amplitude Tf(x, Q 2 ,uj), Eq. ()2.26|) and the gluon 
component of the ?/-meson DA, we have for the gluon part of the form factor 

Q 2 F 9 , . . (Q 2 , uj) = A7[2 ° B ^ j°° due- ut R{u, t) \ C dxx l ~ u x{x - x) 

OOn JO JO 
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+ / dxxx 1 u (x — x) 
Jo 



UJ [(1 + uj)x + (1 — uj)x 
which after some simple calculations becomes 
4tt 2 CB{Q 2 



uj [(1 + uj)x + (1 — uj)x] 



Q F%g*g*(Q ,U) 



3f3 u 



J™ due~ ut R(u, t) |s(4 - u, 2) 2 F 1 (l, 4 - u; 6 - u; - 



(4.21) 



2cj 



+5(4, 2 -u) 2 F 1 (l,4;6-«; T ^) 



2a; 



1 + w 



5(3, 3 -u) 



2 Fi 1 1, 3 - u; 6 - w; — ) + 2 5i f 1, 3; 6 — u; 



1 + w 



1 + w 



1 -co 

1 +CJ 



5(2-^,4)^ 1,2- u; 6- «; 



2cu 
1 + w 



B(2,4-u) 2 F 1 1,2; 6- «; 



2cj 

1 + CJ 



For the ?/-meson on-shell-gluon transition, we find 

47r 2 C5(Q 2 ) /■«> 



.22) 



Q 2 F°, gg *(Q 2 ^ = ±i) = 



/ due- ut R(u, t) 
Jo 



x [5(1, 4 - «) + 5(4, 1 - u) - 5(2, 3 - u) - 5(3, 2 - «)] , 



or equivalently 



Q 2 F° /ggt (Q 2 ,uj = ±l) 



4n 2 CB(Q 2 ) /■«> 



roc / I 

J due~ ut R(u, t) ( — 



u 2 — u 3 — u 



(4.23) 



. (4.24) 



The form factor Q 2 F^ g * g *(Q 2 ,uj = 0) can be calculated by employing the following form 
for the hard-scattering amplitude 



T?(x,Q 2 ,uj)-T$(x,Q 2 ,uj) 



2ir 



a s (Q x) + a s (Q x 



x — x 



1 — uj 2 (x — x) 2 ' 



(4.25) 



which for uj = leads to a very simple expression. Then, it is not difficult to demonstrate 
that 

Q 2 F^ gt (Q\ uj = 0) = f due^R(u, t) [5(4 - u, 2) 



3(3 jo 
— 25(3 — u, 3) + 5(2-m,4)]. 
This latter expression can be recast into the form 

167T 2 C5(Q 2 ^ roc 



(4.26) 



Q 2 5^ w (g 2 ,^ = 0) = 



3,9o 



/ due- ut R(u, t) 
Jo 



2-u 



+ 



3 — u 4 — u 5 — u 



(4.27) 
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The IR renormalon structure of the integrands in Eq. (|4.24jl and (|4.27|) is obvious: they 
have a finite number of IR renormalon poles located at the points uq = 1,2,3 and uq = 
2, 3, 4, 5, respectively. In order to find the IR renormalon structure of the integrand in Eq. 
(|4.22|) . we expand the corresponding hypergeometric functions over (3 = 2u/(l + uS) in the 
region to e (0, 1), providing the following result 



Att 2 CB(0 2 ) r°° 00 



11 1 

k=0 



+B(k + 4-u, 2)] (3 k - - — - [B(k + 2, 4 - u) + B(k + 2 - u, 4)] (3 k 

- [B(k + 3, 3 - u) + B(k + 3 - u, 3)] (3 k+1 } . (4.28) 

Using the identities 

T(A; + 4) 



S(A; + 4,2-m) = 
B(k + 4-u,2 
B(k + 2,4-u) = 
B(k + 2-u,4:) 



(2 - m)(3 - u) . . .(k + 5 - u)' 
1 

~ (fc + 4-u)(fc + 5- u)' 

r(fc + 2) 

(4 -w)(5 -«)...(£; + 5 
6 



5(^ + 3,3-^ 
B(ib + 3 - u, 3) = 



(k + 2 - u)(k + 3 - u)(k + 4 - u)(k + 5 - u) 
T(A; + 3) 



(3-u)(4-m)...(A; + 5-m)' 
2 



(k + 3 - u)(ife + 4 - + 5 - u) 

it is easy to conclude that there is an infinite number of IR renormalon poles, being located at 
the points uq = 2, 3, . . . fc+5; uq = k+4, k+5; uq — 4, 5, . . . A;+5; Mo = k+2, k+3, k+4, k+5; 
uq = 3, 4, . . . k + 5 and «o = + 3, k + 4, + 5. 



V. NUMERICAL ANALYSIS 

We begin this section by comparing the results obtained with the RC method with those 
from the IR matching scheme. In Sect. Ill we have noted that from experimental data only 
values of the nonperturbative parameters fi(2 GeV) and /2(2 GeV) have been extracted. 
We also know that the lowest-order moment integral and hence the parameter entering our 
formulas is ^(2 GeV). Therefore, to make the comparison as clear as possible, we should 
choose the input parameters for the r/-virtual gluon transition in such a way as to determine 
the behavior of the form factor solely with / 2 (2 GeV). This can be easily achieved if we set 
for the Ty'-meson DA parameters 
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Q 2 (GeV 2 ) 

FIG. 3: Scaled rj'g transition form factor —Q 2 F v > gg *(Q 2 ,uj = ±1) vs. Q 2 . The solid line is 
computed within the RC method, the dashed one is found using the IR matching scheme with 
/2(2 GeV) ~ 0.535. The dot-dashed curve is obtained in the framework of the IR matching scheme 
using for /2(2 GeV) the experimental value 0.5. 

Under these circumstances, the gluon component of the vertex function vanishes. To remove 
from the analysis the higher-moment integrals f p (Q), p > 2, we consider only the rj'-meson 
on-shell-gluon transition, i.e., the u — ±1 case. Moreover, we neglect the ~ a 2 s order term 
in Eq. ()3.2|) and set in Eq. (|3.8|) R{u, t) = 1 because in Eq. (|3.19|) ot s is used at the level 
of the one-loop order accuracy. After these simplifications, the FF is given by the following 
expression 

1 Fnr 2 C r°° p~ ut rlii 

Q 2 F v , gg *(Q 2 ,uj = ±1) = / %-™ = -4-kCMQ). (5.1) 

fjQ Jo 1 — u 

Results of our computations are shown in Fig. El where in order to distinguish the various 
curves, these are displayed in the region of Q 2 G [1,4] GeV 2 , whereas at higher Q 2 they are 
close to each other. The RC method and the IR matching scheme 5 both lead almost to 
identical predictions in the entire domain 1 GeV 2 < Q 2 < 25 GeV 2 , provided that in the IR 
matching scheme one uses in expression (J3.19)) the value ^(2 GeV) ~ 0.535 found within 
the RC method (J3.23j) . The curve following from the IR matching scheme deviates from 
the prediction obtained with the RC method only for Q 2 < 1.4 GeV 2 . On the other hand, 
the deviation of that curve, calculated using the experimental value ^(2 GeV) ~ 0.5, from 
the result of the RC method is sizeable in the region Q 2 = 1 4- 2 GeV 2 , reaching ~ 30% at 
Q 2 = 1 GeV 2 . A similar behavior was observed in the calculation of the pion electromagnetic 
form factor, carried out within the context of these methods j6j]|. The difference between the 
solid and the dot-dashed lines in Fig. El is considerably reduced when varying the QCD scale 
parameter A = 0.3 GeV or the experimental value of /2(2 GeV) within their corresponding 
uncertainty limits. But we are not going to make decisive conclusions from these rather 



5 Note that in the calculations within the IR matching scheme, Eq. (|3.19|) with N = 4 has been used. 
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FIG. 4: The quark (left panel) and gluon (right panel) components of the transition form factor 
Q 2 Fi gt » (Q 2 , (J) as functions of Q 2 . The solid curves are obtained using the RC method, whereas 
the broken lines are calculated within the standard HSA. The quark component is computed at 
various oj values by employing the asymptotic DA. In the right panel the correspondence between 
the curves and the input parameter B 2 is B 2 = 8 for the curves labelled 1 and B 2 = 4 for those 
labelled 2. 

model-dependent calculations. Our aim here is to check and demonstrate that the RC 
method and the IR matching scheme predict in the considered region 1 GeV 2 < Q 2 < 
25 GeV 2 almost identical results that do not contradict experiment. 

In order to proceed with the computation of the ?/-meson gluon vertex function and 
explore the role played by the rf gluon content in this process, we have to define the allowed 
values of the free input parameters B\ and B 2 at the normalization point /x 2 , = 1 GeV 2 . 
These parameters determine the shape of the DA's of the quark and gluon components of 
the //'-meson and, in general, they have to be extracted from experimental data or computed 
by non-perturbative techniques. The comparison of the //-meson photon electromagnetic 
transition FF F 7? / 7 (Q 2 ) with the CLEO data leads to the conclusion that the TV-meson DA 
must be close to its asymptotic form with a coefficient B 2 ~ 0.1 deduced in [6|. But this 
conclusion was made by neglecting the contribution of the gluon component of the 7/'-meson 
to the 7/7 transition FF. The investigation of the FF F v r y (Q 2 ) was extended and revised in 
|62^ . In fact, in this work the FF F r) / 7 (Q 2 ) was computed within the RC method by taking 
into account contributions arising due to both the quark and gluon components of the rj' 
meson DA. The comparison with the CLEO data demonstrated that allowed values of the 
Gegenbauer coefficients B 2 (l GeV 2 ) and B 2 (l GeV 2 ) are strongly correlated. They were 
extracted in Ref. jH^] and read 

5f (1 GeV 2 ) = 0, Bf (1 GeV 2 ) G [4, 18], 

B q 2 (l GeV 2 ) = 0.05, B g 2 (l GeV 2 ) G [0, 16], (5.2) 

and 

B q 2 {l GeV 2 ) = 0.1, B 9 2 (1 GeV 2 ) G [-2, 14]. (5.3) 
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FIG. 5: The contribution of the quark a) and gluon b) components of the rf meson to the form 
factor Q 2 F 7llgg *(Q 2 ,uj = ±1) vs. Q 2 . All curves are obtained in the context of the RC method. 
The correspondence between all displayed curves and the parameter B 2 i s: -^f = f° r solid curves; 
B 2 =2 for the dot-dot-dashed curve; B 2 = 4 for the dashed lines; B 2 = 6 for the dot-dashed line, 
and B 2 = 8 for the short-dashed curves. In figure a) curves corresponding only to B 2 = 0, 4, and 
8 are shown. 

In the present paper we select values of the parameters -B|(l GeV 2 ) and B 2 {1 GeV 2 ) that 
obey the constraints ()5.2j) and ()5.3|) . 

It is evident that the non-asymptotic terms in the quark and gluon DA's of the ?/-meson 
proportional to A(Q 2 ) and B(Q 2 ), respectively, affect the asymptotic value of the r/'-meson- 
gluon transition form factor. Therefore, before presenting contributions from these terms, 
it is instructive to study the asymptotic FF itself. In the left panel of Fig. HJ we depict 
the ?/-meson - virtual gluon transition FF as a function of the gluon virtuality Q 2 . For the 
asymptotic DA the quark component of the form factor coincides with the full one. In the 
same figure the predictions obtained within the standard HSA are also shown. One sees 
that in the domain 1 GeV 2 < Q 2 < 25 GeV 2 the standard pQCD results get enhanced by 
approximately a factor of two due to power corrections. A similar conclusion is valid also 
for the gluon component of the form factor (right panel in Fig. 01 computed by employing 
the r/-meson DA's). 

Here some comments concerning the accuracy of the numerical computations are in order. 
In Fig. |3] (left panel) and in the following ones, the curves obtained within the RC method, as 
a rule, require a summation of an infinite series. In real numerical computations we truncate 
such a series at some k = K max . Naturally, the question arises about the convergence rate of 
this series. Let us explain this problem by considering Fig. 0] as an example. The solid curve 
with uo = 0.5 in Fig. 0] (left panel) has been found within the RC method by employing Eq. 
(I4.12JI . This expression contains a series with factorially growing coefficients. For definiteness 
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FIG. 6: The form factor —Q 2 F >n i gg *(Q 2 ,uj = ±1) computed using the RC method a) and the 
standard HSA b). The correspondence between displayed curves and the parameter 7?f is: = 
for solid curves; -Bf = 2 for the dot-dot-dashed curves; -Bf = 4 for the dashed lines; -Bf = 6 for the 
dot-dashed lines. 

we analyze the term 

/*oo °° 

/ due- ut R(u,t)J2B(k + 2,2-u)p k 

£ r due^R(u,t)- w/V 2 )*.-^ S P k - (5-4) 

f— Wo (2 — u)(6 — u) . . . (k + 3 — u) 



k=0 



The expansion parameter (3 = 2u/(l + uj) in Eq. (|5.4|) at the point u; = 0.5 is equal to 
(3 = 2/3. Below, we write down the values of the Gamma function T(k + 2) = (k + 1)! and 
also those of the product of (5 k with the principal value of the integral 

poo 

Hk)= due- ut R(u,t)- £ (5.5) 

V ' Jo K ' J (2-u)(3-u) ...{k + 3-u) K ' 

for k = 0, 5 and 10 

r(2) = i, r(7) = 720, r(i2) = 3.99168 • 10 7 , 

and 

1(0) ~ 0.09041, 7(5) ~ 2.44061 ■ 10" 6 , 7(10) ~ 1.31568 ■ 10" 12 . 
As a result, the corresponding terms in the sum given by Eq. (J5.4)) take the values 

0.09041, 1.75724 • 10~ 3 , 5.25177 ■ 10~ 5 , 
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FIG. 7: The quark a) and gluon b) components of the FF at fixed total gluon virtuality 
Q 2 = 4 GeV 2 as functions of the asymmetry parameter u. The RC method is employed. The 
correspondence between depicted lines and the parameter is the same as in Fig. |21 

respectively. The calculations above have been performed at Q 2 = 2 GeV 2 . At the momen- 
tum transfer Q 2 = 20 GeV 2 we get 

0.03898, 7.59618- 10" 4 , 4.42997 ■ 10~ 5 . 

One observes that the convergence rates of the numerical series are high and that we can 
therefore truncate them, rule, at K max = 20. 

We have analyzed the impact of the various DA's of the f/-meson on the rj'g transition 
form factor. The quark component of the FF is stable for different values of i?| £ [0,8]. 
It is difficult to distinguish the corresponding curves and therefore in Fig. 5a we can only 
plot some of them. In contrast, the gluon component of the form factor demonstrates a 
rapid growth with i?| (Fig. 5b). As a result, due to different signs of the quark and gluon 
components of the space-like vertex function, the total vertex function Q 2 F v ' gg *(Q 2 ,u = ±1) 
for 7^ lies below the asymptotic one (Fig. 6a). Comparing the predictions derived within 
the RC method with those following from the standard HSA (Fig. 6b), we see a quantitative 
difference between corresponding curves. 

The dependence of the quark and gluon components of the form factor on the asymmetry 
parameter uj at fixed Q 2 and for various DA's of the //-meson are shown in Fig. [7| Because 
F 1] /g*g*(Q 2 , uj) is symmetric under the exchange to <-> —uj, we present our results in the region 
< uj < 1 only. We have just demonstrated in Fig. El that the quark component of the FF 
F^, gg *(Q 2 ,uj = ±1) is not sensitive to the use of various ?/-meson DA's, employed in our 
calculations. This is valid also for its behavior as a function of uj (Fig. 7a). In accordance 
with our computations, an effect of the chosen parameter B\ on the gluon component of the 
FF is considerable in the whole range of uj G [0, 1]. 

The magnitude of the quark and gluon components of the form factor for a given DA 
depends on the total gluon virtuality Q 2 (Fig. |HJ). In this case both the quark and gluon 
contributions to the FF demonstrate sensitivity to the fixed value of Q 2 . 
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FIG. 8: The quark a) and gluon b) components of the form factor Q 2 F„i g * g * (Q 2 , to) at fixed Q 2 
vs. oj. The RC method is used. The momentum scale for the solid curves Q 2 = 10 GeV 2 , for the 
dashed ones Q 2 = 6 GeV 2 , and for the dot-dashed curves Q 2 = 4 GeV 2 . 

The features of the quark and gluon components of the r]'g* transition FF described above 
determine the behavior of their sum as a function of the asymmetry parameter u and, as a 
result, we get the picture shown in Fig. 9a. Owing to the gluon component, the form factor 
F n i g * g *(Q 2 ,uj) depends on the 7/-meson DA used in the calculations. For comparison, in Fig. 
9b the curves found within the standard HSA are also shown. An enhancement of about a 
factor of two of the vertex function due to power corrections is evident. 

The 7]'g* transition FF as function of the first gluon virtuality Q 2 at various fixed values 
of the second one, Q 2 ,, and for different DA's is plotted in Fig. [TU1 

As we have noted in Sect. Ill, the principal value prescription, adopted in this work to 
regularize divergent integrals, produces higher-twist ambiguities. It is important to clarify to 
what extent these ambiguities may alter our predictions. To illustrate this effect, we depict 
in Fig. ITTl as an example, the scaled rj'g transition FF Q 2 F v / gg * (Q 2 , u = ±1). The higher- 
twist ambiguities lead approximately to ±15% shifts relative to the previous RC result (line 
2). On the other hand, comparing the standard pQCD prediction (solid line 1) with the 
lower dashed line (that denotes the RC prediction with a negative higher-twist ambiguity), 
we see that still the power corrections provide an enhancement of the standard pQCD result 
by a factor ~ 2 in the region Q 2 ~ 1 - 2 GeV 2 and by a factor ~ 1.6 at Q 2 = 25 GeV 2 . This 
latter effect is connected not only to larger uncertainties, but also reflects the general trend 
of the FF to reach its asymptotic value, i. e., the line 1 in the limit Q 2 — ► oo. 

VI. CONCLUDING REMARKS 

In this paper we have evaluated power-suppressed corrections ~ 1/Q 2n , n — 1, 2... to 
the space-like //-meson virtual-gluon transition form factor Q 2 F v / g * g *(Q 2 ,u). To this end, 
we have employed the standard hard-scattering approach and the running-coupling method 
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FIG. 9: The FF —Q 2 F v i g * g *{Q 2 ,uj) obtained by employing the RC method a) and the standard 
HSA b) as a function of uj. The expansion coefficients are 5| = for the solid lines, -Bf = 4 f° r 
the dashed lines, and -Bf = 8 for the dot-dashed lines. 

in conjunction with the infrared-renormalon calculus. In the calculations, both the quark 
and the gluon distribution amplitudes of the rf meson have been taken into account. In 
these model DA's only the first non-asymptotic terms have been retained and the values of 
the input coefficients -Bf(l GeV 2 ) and £?f (1 GeV 2 ) extracted from the analysis of the CLEO 
data on the 7/7 electromagnetic transition FF have been employed. 

In order to apply the RC method to the considered process, the hard-scattering ampli- 
tudes of the corresponding subprocesses have been generalized in such a way as to preserve 
the symmetry properties of both the hard-scattering amplitudes and the transition form fac- 
tor itself under the replacements x <-» x and uj <-> — uj. In the computations within the RC 
method, the Laplace transformed expression for the running coupling has been employed. 
The Borel resummed form factors, obtained this way, have been regularized by means of the 
principal value prescription. Various limits of the general expression Q 2 F r] ' g * g *(Q 2 ,uj) have 
been found. 

Our expressions for the vertex function F^ g * g * (Q 2 , uj) , found within the standard HSA, 
are in agreement with corresponding predictions made in Ref. j26j (up to a conventional 
sign factor). It has been demonstrated that (see Eq. (|2.27Jl ) both the quark and the gluon 
components of the vertex function are symmetric under the exchange uj <-> —uj, as they must 
be owing to the Bose symmetry of the final gluons. The RC method has provided us a tool 
for estimating power corrections to the vertex function F v i g * g *(Q 2 ,u). These corrections are 
implicitly contained in the pQCD factorization formulas (j2.14|) . (j2.15p and originate from the 
end-point x — > 0; 1 regions. It is clear that such corrections cannot been taken into account 
in the standard pQCD approach by freezing the renormalization scale fi 2 R and ignoring its 
dependence on the longitudinal momentum fraction x. As an important consistency check, 
we have proven that the results obtained with the RC method in the asymptotic limit 
Q 2 — * 00 reproduce the standard pQCD predictions for the vertex function. This provides 
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Q, 2 (GeV 2 ) Q, 2 (GeV 2 ) 

FIG. 10: The form factor —Q 2 F rj i g * g * (Q1,lo) for two different DAs of the rj meson and at different 
fixed values of Q 2 ,', viz., for the solid curves Ql = 1 GeV 2 , for the dashed curves Ql = 5 GeV 2 , for 
the dot-dashed curves Q\ = 10 GeV 2 , and for the dot-dot-dashed curves Q 2 = 25 GeV 2 . 

further justification for the treatment of the hard-scattering amplitudes ()2.29|) . (|2.30j) . and 
the symmetrization procedure employed in the running coupling method. 

The presented numerical analysis shows that power corrections considerably enhance the 
standard pQCD predictions for the form factor in the explored region 1 GeV 2 < Q 2 < 
25 GeV 2 , though other sources, not considered here, may also give rise to power corrections. 
Our investigations demonstrate that the quark component of the form factor at fixed B\ = 
0.1 is practically stable for several values of B\ = 0, 2, 4, 6, 8. Contrary to this, the 
gluon component of the FF is sensitive to the adopted value of B\. As a consequence, the 
r]'g* transition FF was found to depend on the gluonic content of the rj' meson. In the 
considered region the gluon contribution reduces the absolute value of the space-like form 
factor Q 2 F ri ,g* g *(Q 2 ,uj). 

As is true of any technique for calculating power corrections, the theoretical framework 
elaborated in the present work makes assumptions about the regularization of endpoint di- 
vergences. There are of course other possibilities. Nevertheless, we believe that our method 
is useful in pQCD analyses of the B meson exclusive decays and heavy-to-light B — > ir, p 
transition form factors in the domain of moderate momentum transfers. A generalization of 
the RC method to describe time-like transitions as well as its combination with resumma- 
tion techniques to include Sudakov logarithms will be the subject of separate investigations. 
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APPENDIX A 



In this Appendix we collect expressions for the //-meson quark and gluon DA's for rif 
The parameters (|2.38j) which determine these DA's assume the following values 



2 35 2 2 48 2 107 

'99 3 ' '99 **■)!+— g ) 1- — g ) 

19 1 

q xiJ q 1 
P2 ~ T> ^2 



5 ' ' z 102 

It is evident that the elements of the anomalous dimensitions matrix 7^ and 7^ are rif 
independent and have the values shown in Eq. (|2.38|) . 

In the rif — 4 case the DA's 4> q (x, Q 2 ) and <j) 9 (x, Q 2 ) of the 7/-meson have the form ()2.39j) 
as well, the only difference being in the functions A(Q 2 ) and B(Q 2 ), which now are defined 
by the expressions 



A{Q 2 






V(Q 2 )\ 




V«s(/^o)/ 


48 

) 75 + 5Bg 


/« s (g 2 )\ 



107 
7S 



107 
75 



5(Q 2 ) = 1951 +55fP^ . (A.l) 



In the previous sections all results have been written down for rif = 3 valid for a mo- 
mentum transfers 1 GeV 2 < Q 2 < 2 GeV 2 . For momentum transfers in the range 
2 GeV 2 < Q 2 < 25 GeV 2 , the choice rif = 4 has to be employed. In this case, the 
expressions that determine the quark component of the transition FF remain unchanged, 
except for the function A{Q 2 ) (and /5o, R(u,t)), which should be taken from Eq. (jA.l|) . The 
expressions for the gluon component of the FF have to be rescaled by a factor 3/4 — apart 
from the replacement of the function B{Q 2 ) (and /3 , R(u,t)) — since the hard-scattering 
amplitudes T^Jx, Q 2 , uj) explicitly depend on rif (|2.30|) . 
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